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Abstract 

Overshooting occurs in stars when convective elements penetrate into adjacent 
radiative zones. In the Sun, it leads to the so-called 'tachocline' at the base 
of the outer convection zone and this region is becoming a key ingredient of 
the standard solar dynamo model as strong toroidal magnetic fields may be 
generated there. However this overshoot is not predicted by the mixing-length 
II , theory of convection where convective elements must stop at the border of a 
O ' convectively unstable region. I will review the main properties of this convective 
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overshoot in stellar interiors, with in particular its subtle dependence with the 
thermal diffusion, and will present the most recent results obtained from 2-D 
and 3-D direct numerical simulations (DNS) of penetrative convection. 



Session: Physics and uncertainties in massive stars on the MS and close to it 

m : 

"sf , Introduction 

OO ■ 

, In the standard model of thermal convection based on the mixing-length phe- 
t-H ■ nomenology (Bohm-Vitense 1958), the penetration of convective elements into 
[ surrounding stable layers is not possible. Indeed, the mixing-length velocity v 
of a convective blob is related to its temperature contrast ST by 



jrrp 

v 2 oc —gl with <STocV- V ad , (1) 

where g and I denote the gravity and the mixing length, respectively, and V 
C3 | and V a d the usual temperature 'nabla'. Following Schwarzschild's criterion, 
the interface that separates a convective zone from a radiative one corresponds 
to V = W a d therefore the blob must stop here. However, the mixing-length 
formalism is based on an acceleration budget and setting the blob acceleration 
to zero does not mean that its velocity is itself zero, i.e., the mixing-length 
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theory neglects inertia. In fact, the blob penetrates into the stable layer over a 
small distance, then it is slowed down by the buoyancy braking and finally stops. 
This small penetration is not anecdotal because it has important consequences 
in stellar physics: 

- it leads to different evolutionary tracks in the HR diagram as massive stars 
live longer on the main sequence due to the additional mixing induced by the 
overshooting above the convective core (Rosvick & Vandenberg 1998, Perryman 
et al. 1998). 

- in the standard a — f2 model of the solar dynamo, the regeneration of poloidal 
magnetic fields into toroidal ones takes place in the tachocline where vertical 
and latitudinal differential rotation is expected to be most efficient (Parker 
1955, Brandenburg & Subramanian 2005). 

- heliosismology showed that the solar core rotation is almost rigid (Brown et 
al. 1989). Internal gravity waves propagating in stellar radiative zones are 
suspected to play the main role in this angular momentum redistribution (Zahn 
et al. 1997, Talon & Charbonnel 2003). The most wide-spread excitation 
model involves penetrative convection from neighboring convection zones as 
strong downdrafts extend a substantial distance into the adjacent stable zones 
so that internal gravity waves can be randomly generated (Dintrans et al. 2005). 

Modelling of overshooting 

Assuming that the temperature gradient in the overshoot layer is close to the 
adiabatic one of the next convection zone, Roxburgh (1978) investigated the 
overshoot of convective cores and derived what it is now called the 'Roxburgh's 
integral constraint': 

£(L-L rad )d(J^)>0, (2) 

where the radius r c corresponds to the edge of the overshoot region, while L 
and L ra d denote the total and radiative luminosity, respectively. The overshoot 
layer corresponds in this relation to the region where L < L ra d (i.e., negative 
contribution) while the convection zone means L > i rac j (i.e., positive contri- 
bution). However, this relation gives an upper limit of the overshoot extent 
because it does not take into account the dissipation acting on the penetrating 
elements. Indeed, the viscous dissipation T> v enters in the heat equation and 
affects the blob dynamics such that the correct relation reads now (Roxburgh 
1989): 

£(L - L M )4 (I) = [ 4^5* with V. ~ , (£)* . P) 
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where v denotes the kinematic viscosity. As the viscous dissipation is always 
positive, this term acts to decrease the negative part of the integrand in the 
left-hand side and thus the radius r c . 

Following Roxburgh's investigation, Zahn (1991) showed that the penetra- 
tion extent also depends on the value of the Peclet number Pe associated with 
the convective elements, that is, 



where \ denotes the thermal diffusivity. Small Peclet's numbers mean that 
the radiative diffusion is stronger than the advection and the convective blob 
rapidly loses its identity compared to the surrounding medium. As a conse- 
quence, it does not feel any more the buoyancy braking and travels over long 
distances in the radiative zone (Zahn called this regime 'overshoot'). On the 
other hand, a blob with a large Peclet number does not thermalize rapidly with 
its surrounding (i.e., no radiative losses) and keep its identity: it is then strongly 
slowed down by the buoyant force and its penetration is weaker (the so-called 
'penetration' regime). Penetration leads to an adiabatic stratification below the 
convection zone because of the induced entropy mixing, whereas overshoot has 
no significant influence on the local stratification. 

Assuming that the velocity W of the penetrating motions obeys the usual 
mixing-length scaling of thermal convection, Zahn derived the following rela- 
tions for both the subadiabatic extent of a convective region and the width of 
the thermal boundary layer: 



where / denotes a filling factor of convective plumes (i.e., fraction of the area 
occupied by the plumes) and H p the pressure scale-height at the bottom of 
the convection zone. Applied to the Sun, these two scalings lead to a pene- 
tration extent of about 0.2 - 0.3_ff p (~ 15 000 km) and an overshoot layer 
~ 1 km. These scalings seem however to overestimate the amount of con- 
vective penetration inferred by helioseismology as the study of the oscillations 
in the asymptotic laws of solar acoustic modes rather predicts an upper limit 
of ~ O.lHp for the extent of the whole convective overshooting (Roxburgh & 
Vorontsov 1994; Christensen-Dalsgaard et al. 1995; see also Aerts et al. 2003 
for the case of a massive star). 
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Numerical simulations of overshooting 
The bad news 

The rapid development of supercomputers has lead several groups to under- 
take direct numerical simulations of convective overshooting. However, many 
numerical difficulties appear when one tries to simulate stellar convection: 

- the molecular viscosity in stellar interiors is small, typically v ~ 1 cm 2 /s 
(Parker 1979). Granulation cells in the upper convection zone of the Sun 
correspond to a typical size and velocity of about L = 1000 km and U = 1 
km/s, respectively. It leads to Reynolds' number of about Re = UL/v ~ 10 11 , 
meaning that solar convection is a highly turbulent and nonlinear phenomenon. 

- we know from Kolmogorov's theory of turbulence (Kolmogorov 1941) that 
the ratio between the largest scale (i.e., the injection scale L) and the smallest 
one (i.e., the dissipation scale id) behaves as L/td ~ -Re 3 / 4 such that the 
solar dissipation scale is about one centimeter. In other words, one needs 
(at least) 10 8 gridpoints in each spatial direction to well reproduce both the 
injection of energy in the largest structures and its viscous dissipation in the 
smallest ones (the turbulent cascade scenario). As this number of gridpoints is 
clearly beyond the scope of today supercomputers, several numerical techniques 
have been developed to overcome this restriction: (i) instead of doing Direct 
Numerical Simulations (DNS), one can perform Large-Eddy Simulations (LES) 
where the small-scale dynamics is modelled using a sub-grid scale approach 
(Lesieur 2008); (ii) one can also artificially increase the dissipation scale by 
adding a large-scale viscosity in the whole computational domain, the price to 
pay being that motions are essentially laminar [Re ~ 10 2 • • • 10 4 ). 

- finally, another problem concerns the different time scales that coexist in a 
convection zone. Indeed, the local turnover time scale of convection is roughly 
given by: 

t ~ H. p IU R ms, (6) 

where Urms is the typical (turbulent) RMS velocity of convective eddies. Ap- 
plied to the solar convection zone, it leads to turnover time scales of about the 
day at the surface to several months at the bottom (Spruit 1974). As a con- 
sequence, simulations should reproduce motions whose dynamical time scales 
can vary by two orders of magnitude in the computational domain. Moreover, 
the thermal time scale is also much more longer than these turnover time scales 
by orders of magnitude as, e.g., the solar thermal flux expressed in dynamical 
units is very small at the bottom of the convection zone (F /(/9c 3 ) ~ 10~ n ). 
This ratio is in essence equal to the ratio of the dynamical to thermal time 
scales, meaning that realistic simulations should be integrated on a very long 
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Figure 1: Penetration extent A v.s. stability parameter S. Left: in the 2-D case 
where both the penetration and overshoot regimes exist (from Hurlburt et al. 1994). 
Right: in the 3-D case where only the overshoot regime remains (from Brummell et 
al. 2002). 



time to make sure that a thermal-relaxed state is reached. It is not possible in 
practice to do that with a decent spatial resolution and the thermal fluxes are 
commonly overestimated by orders of magnitude in numerical simulations. 



DNS of convective overshooting in local cartesian boxes 

Following Hurlburt et al. (1984), the majority of DNS done so far use polytropic 
solutions for the internal structure of the star. Indeed, polytropes with p oc 
T m (m being the polytropic index) allow to easily specify the strength of the 
convective instability as the hydrostatic equilibrium implies in that case 

d S = i ( 7) 

where g is the downward-directed gravity, z the altitude and i?„ the perfect gas 
constant. For a monatomic gas with an adiabatic index 7 = 5/3, this relation 
shows that the temperature gradient is equal to the adiabatic one (= —g/c p ) 
when the polytropic index is m ac i = 1/(7—1) = 3/2. As a consequence, 
a stable layer corresponds to to > m ac i and a convectively-unstable one to 
m < m ad . 

Hurlburt et al. (1986, 1994) used this model to study the overshooting in 
2-D simulations of compressible convection for a solar case, that is, a radiative 
zone with polytropic index 777,3 > 3/2 located above a convective zone with 
polytropic index 7772 < 3/2. In particular, they studied the influence of the 
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relative stability parameter S on the penetration extent A, where S = — (m a d — 
m3)/(m a d — TO2), and found the following scalings (Fig. 1, left): 

A oc S^ 1 for small S-values: penetration regime, 

(8) 

A oc S 1//4 for large S-values: overshoot regime. 

However, these scaling laws have not been fully confirmed in the 3-D case by 
Brummell et al. (2002) as only the overshoot regime A oc is recovered 

(Fig.l, right). The reason simply lies in the smaller filling factor / for the con- 
vective plumes compared to the 2-D case (Eq. 5): the spaced downdrafts that 
penetrate into the stably stratified layer are not enough to maintain an adia- 
batic stratification below the convection zone and only the thermal boundary 
layer exists. 




Figure 2: Vorticity isocontours in a 2-D star-in-a-box simulation of penetrative con- 
vection for a solar-like star (Dintrans & Brandenburg 2009). 



Another interesting 2-D and 3-D simulations of penetrative convection are 
the ones done by Roxburgh & Simmons (1993) and Singh et al. (1998). Indeed, 
a more realistic model of the star interior has been considered in the 2-D case 
through a temperature-dependent radiative conductivity. The main result is 
that the viscous dissipation contribution in the RHS of Roxburgh's integral 
constraint (3) decreases with the Prandtl number. It suggests that in the 
astrophysical limit of very small Prandtl numbers (Pr ~ 10~ 6 ), the stable and 
unstable contributions in the LHS of Roxburgh's integral strictly balance each 
other and are sufficient to obtain a good estimate for the overshoot extent 
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(Roxburgh & Simmons 1993). In the 3-D case (Singh et al. 1998), a sub- 
grid scale modelling of the small scales led to the confirmation of the scaling 
relationships between this extent and both the imposed bottom flux and the 
vertical velocity of convective downd rafts (see Eq. 5). 

Global DNS 

In the last decade, global direct simulations of overshooting have greatly im- 
proved following two numerical approaches: 

- the spherical shape of the star is reproduced using a pseudo-spectral code 
where physical fields are projected onto spherical harmonics. Browning et al. 
(2004) used this method with the 3-D anelastic ASH code to study the core 
overshooting in a massive A-type star of 2M Q , where only the inner 30% by 
radius of the star has been considered (convective core and some of its surround- 
ing radiative envelope). Although the achievement of the thermal relaxation 
in these simulations remains questioning, they found that the convective core 
adopts a prolate shape while overshooting is larger at the equator than at the 
poles, yielding an overall spherical shape to the whole central region. 

- the star is embedded in a topological^ Cartesian domain and the sphericity 
is reproduced through radial-dependent forcing or damping profiles. This novel 
approach has been developed by Dintrans & Brandenburg (2009) and first 
results in 2-D are promising as both the penetration and overshoot regimes are 
observed, whereas internal gravity waves are efficiently excited in the central 
radiative zone (Fig. 2). 

Conclusion 

Convective overshooting is today included in stellar evolution codes and evo- 
lutionary tracks of massive stars reproduce satisfactorily the observations of 
clusters when an overshooting above convective cores is assumed. Theoreti- 
cal studies showed that the overshoot region is composed of two nested layers 
whose the size depends on the value of the Peclet number of penetrating mo- 
tions. Numerical studies in 2-D and 3-D of penetrative convection confirmed 
this picture, except that only the overshoot layer (i.e., the thermal boundary 
one) exists in 3-D due to the weak filling factor of convective plumes. These 
local and global direct numerical simulations are however still quite far from 
realistic values of stellar interiors. Flows are less turbulent than in reality and 
the overshoot regime dominates. Further progress are clearly needed in the 
numerical side of this problem (e.g., development of large-eddy simulations), in 
its theoretical modelling (e.g., turbulent closure models; Kupka & Montgomery 
2002) and in the observational signatures of overshooting from asteroseismol- 

ogy- 
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